High Magnetic Shear Gain in a Liquid Sodium Stable Couette Flow Experiment 

A Prelude to ru a — Q Dynamo 



Stirling A. Colgate^, ^, Howard Beckley^, (deceased), Jiahe Si^, Joe Martinic^, David Westpfahl^, James 
Slutz^, Cebastian Westrom^, Brianna Klein^, Paul Schendel^, Cletus Scharle^, Travis McKinney^, Rocky 
Ginanni^, Ian Bentley^, Timothy Mickey^, Regnar Ferrel^, and Hui Li^ Vladimir Pariev^, John Finn^ 
^ T-2, MS B-227, Los Alamos National Laboratory, Los Alamos, NM 87545; colgate@lanl.gov 
^Department of Physics, New Mexico Institute of Mining and Technology, Socorro, NM 87801 and 
^T-5, Los Alamos National Laboratory, Los Alamos, NM 87545; 

The Q-phase of the liquid sodium a-Q dynamo experiment at NMIMT in cooperation with LANL 
has been successfully demonstrated to produce a high toroidal field that is ^ 8 x , where Br is 
the radial component of an applied poloidal magnetic field. This enhanced toroidal field is produced 
by the rotational shear in stable Couette flow within liquid sodium at magnetic Reynolds number 
Rm ^ 120. The small turbulence in stable Taylor-Couette flow is caused by Ekman flow with an 
estimated turbulence energy fraction of {6v/vf - 10"^ This high Q- gain in low turbulence flow 
contrasts with a smaller Q-gain in higher turbulence shear flows. This result supports the ansatz 
that large scale astrophysical magnetic fields can be created by semi-coherent large scale motions in 
which turbulence plays only a smaller diffusive role that enables magnetic flux linkage. 
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Major efforts are spent world-wide on astrophysical 
phenomena that depend upon magnetic fields, e.g., plan- 
etary, solar, and stellar magnetic fields. X-rays, cosmic 
rays, TeV gamma rays, jets and radio lobes powered by 
active galactic nuclei (AGNs), and in the interpretation 
of Faraday rotation maps of galaxy clusters. Yet, so far, 
there is no universally accepted explanation for the in- 
ferred lar age-scale magnetic field. A process, called the 
a — Q dynamo [10, 16, 23 , has been proposed which 
involves two orthogonal conducting fluid motions, shear 
and helicity. When the two motions are comparable, it 
is often described as a stretch, twist, and fold or "fast" 
dynamo, and supposedly can be produced by turbulence 
alone [18], [33]. The problem is that a turbulent dynamo 
must create these two orthogonal motions turbulent mo- 
tions. Fluid turbulence maximizes entropy and the diffu- 
sion of the flux of momentum [31 , so that a field twisted 
one way by turbulent eddies at one moment of time may 
be partially twisted the opposite way with only a net 
bias, and hence partially diffusive. 

One possible way to achieve dynamo is to use naturally 
occurring large scale near-stable rotational shear flows 
(as in AGN accretion disks and in stars) in combination 
with transverse, transient, rotationally coherent plumes. 
Such plumes may be driven either by star-disk collisions 
in AGNs [21] or large scale (density scale height) convec- 
tive elements in the base of the convective zone in stars 
[31 US] [32] . These flows can stretch and wind-up an em- 
bedded, transverse magnetic flux through a large number 
of turns (the Q effect). In both types of shear flows the 
turbulence is expected to be relatively small because of 
the stability imposed by either an angular momentum 
gradient or an entropy gradient. (Positive work must be 
done to turn over an eddy in such a gradient.) The ad- 
vantage of a large f^-gain is that the corresponding a-gain 
(the helicity deformation necessary to convert a fraction 



of the toroidal flux back into poloidal flux) can be corre- 
spondingly smaller. Parker, the originator of the a — Q 
dynamo concept for astrophysical dynamos, invoked cy- 
clonic motions to produce the helicity, [24]. Moffatt jl6] 
countered that a typical atmospheric cyclone makes very 
many turns before dissipating, therefore making the sign 
of the resulting poloidal flux incoherent. Moffat sought 
motions where the rotation angle was finite, ~ 7r/4. Wil- 
lette [32 and the first author suggested that buoyant 
plumes were a unique solution to the problem of coher- 
ent dynamo helicity. Beckley et al. [1 demonstrated 
this property experimentally of finite, ~ 7r/4, coherent 
rotation of a driven plume before it dissipated through 
turbulent diffusion in the background fluid. 

Two liquid sodium dynamo experiments have pro- 
duced positive exponential gain, but the flows were con- 
strained by rigid walls unlike astrophysical flows [HI [30] . 
The rigid walls restrict turbulent eddy size by the dis- 
tance from the wall, log-law of the walls [T?, thereby 
producing low turbulence. The three recent experiments 
use the Dudley- James [6 or Von Karman flow (counter 
rotating flow converging at the mid-plane and driven by 
two counter rotating propellers or two vaned turbine im- 
pellers, respectively, [Zl |13l |Il[20l|25l |26l [29] ) . Turbu- 
lence is induced by the Helmholtz instability at the shear- 
ing mid-plane. This combination of coherent motions and 
unconstrained shear-driven (relatively high level) turbu- 
lence resulted in a maximum ^l-gain of only ~ x2. Rec- 
ognizing the enhanced resistivity of the mid-plane tur- 
bulence, the team at Wisconsin added a mid-plane baffle 
to reduce the turbulence, following which the l]-gain in- 
creased to x4 [27 . The von Karman Sodium 2 (VKS2) 
experiment in the same geometry produced exponential 
gain [17 ; However, the dynamo action was explained not 
by turbulence but primarily by the production of helic- 
ity by the large coherent vortices produced by the radial 
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rigid vanes of the impeller [TT] . (Ferro magnetism added 
additional gain to this source of helicity.) 

The New Mexico a — Q dynamo (NMD) experiment 
is a collaboration between the New Mexico Institute of 
Mining and Technology and Los Alamos National Lab- 
oratory m [5] . It is designed to explore the possibility 
in the laboratory using coherent fluid motions in low- 
turbulence high shear Couette flow in the annular vol- 
ume between two coaxial cylinders, R2/R1 = 2 rotating 
at different angular velocities (see Figure 1), closely anal- 
ogous to natural fluid motions that occur in astrophysical 
bodies [2j|21]|22]. The results of the first phase ((^-phase) 
of this two-phase experiment are presented here. 

Fig. 1 also indicates the Ekman flow, a thin fluid layer 
flowing along the boundaries of the annular volume and 
including a small radial return flow from inner to outer 
cylinder through the fluid volume. The Ekman flow pro- 
duces both a torque and a small but finite level of turbu- 
lence [9 . This turbulence adds a small turbulent resis- 
tivity to the resistivity of metallic sodium. Five pressure 
sensors are distributed radially at one end. They mea- 
sure the pressure distributions, giving information of the 
actual rotational flow angular velocity distribution. 
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FIG. 1. A schematic drawing of the fluid flow of the exper- 
iment shows the inner cylinder of Ri — 15.2 cm rotating at 
Qi/27r = 68Hz relative to the outer or confining cylinder of 
R2 = 30.5 cm at Q2 = Qi/A. Liquid sodium in the annular 
space undergoes Couette rotational flow. In addition, a thin 
layer in contact with the cylinder walls and end walls under- 
goes Ekman flow, circulating orthogonally around and then 
through the Couette flow. Pressure sensors in the end wall 
and the magnetic probe housing are also shown schematically. 

For the magnetic measurements, a 50 kW AC induc- 
tion motor (20kW is used for stable Couette flow) drives 
a gear train with clutches and power take-off that ro- 
tates the two cylinders at the fixed ratio of Qi = 41^2 • 
The outer cylinder can also be disengaged from the gear 
train by a clutch, and a DC motor is used to accelerate or 
brake the outer cylinder independently from the driven 
inner cylinder. This allows different Q1/Q2 ratios to be 
explored. The DC motor housing (stator) is mounted 
on bearings. A torque arm with two force sensors (H — ) 



connects the motor stator to ground, so that the torque 
on the DC motor can be measured separately from the 
drive of the inner cylinder. This arrangement allows us 
to measure the torque between the two cylinders due to 
the Ekman flow along the surface of the end plates which 
rotate more slowly at 1^2- 

In particular, when the inner cylinder is driven at 
higher speed by the AC motor, the Ekman fluid torque 
tries to spin up or accelerate the outer cylinder. Two 
torques counteract this acceleration: 1) the friction in 
the bearings that support the rotation of the outer cylin- 
der and 2) the torque on the DC motor when used as 
a generator, or brake. (The generated power is dissi- 
pated in a resistor.) In Fig. 2 (left) the crosses are the 
measurements of the braking force exerted by the DC 
motor torque arm; the dashed line is the calibrated bear- 
ing torque (measured by disengaging the inner cylinder 
drive and rotating the outer cylinder with the DC mo- 
tor alone). The sum of these two torques is equal to the 
Ekman fluid torque spinning up the outer cylinder. 
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FIG. 2. (left) The lower crosses are the measured DC genera- 
tor torques; the dashed line is the measured bearings torque, 
and the sum of the two is the top solid curve, the torque 
due to the Ekman flow, (right) Shown is the measured and 
theoretical pressure distribution. The pressure measurement 
starts at Ri and extends to R2. The finite pressure at Ri, 
Pi :^ 100 psi is due to the centrifugal pressure of rotating 
fluid in the narrow gap between the inner cylinder rotating 
at Qi and the end wall rotating at Q2 (dot-dash curve). This 
thin end-gap is filled with rotating fluid, which extends from 
the shaft seal (where P = 1 atm) at radius, R3 ^ Ri/3 to 
Ri at which radius the Couette flow in the annular volume is 
initiated. From there, the pressure in the annular volume is 
extrapolated either as Q oc R~^ or R~^^^ . The measurements 
fall along the R~^^'^ pressure distribution curve indicating less 
shear than the ideal Couette flow. 

The experiment was originally designed upon the sup- 
position that the primary torque between the two cylin- 
ders would be due to the radial Ekman flow at the end 
walls and would be small, ~ 1/10, compared to pipe flow 
wall friction stress, r = pCnv'^' In Fig. 2 (left) the mini- 
mal Ekman flow torque occurs at (^1/(^2 4, somewhat 
less than the limit of stable Couette flow. The torque 
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value is 2 X 10^ dyne-cm, close to the approximation 
that each Ekman layer thickness [9^ is A :^ RxjRe^l^ :^ 
5 X 10~^ cm and that the average radial flow velocity 
within the layer is Vr — (yt\R\)l2. In Fig 2 (right), 
the measured pressures are compared to the calculated 
pressures corresponding to two different angular velocity 
power laws. The upper (dashed) curve corresponds to 
ideal, maximum shear, stable Couette flow, or ^7 cx r~^, 
but the experimental points follow the lower solid curve 
corresponding to ex r~^/^. The Ekman flow torque 
has distorted the angular velocity profile and reduced the 
shear relative to the ideal Couette flow. 

Turbulent flow at such high Reynolds number Re 
10^ is still well beyond current simulation capability. The 
Ekman flow is a flux of fluid of reduced angular momen- 
tum deposited at the inner cylinder. The torque reducing 
this angular momentum is friction with the end walls, 
which in turn reduces the angular velocity of the Cou- 
ette flow in the annular volume. Unstable flow, turbulent 
friction with the inner cylinder surface, counteracts this 
torque by speeding up this flow. The torque in the Cou- 
ette volume is a constant, independent of radius and axial 
position. Then this shear stress is maintained constant 
in two laminar sub layers [28 with Re 100 and in two 
turbulent boundary layers in a log-law-of-the-walls solu- 
tion [l^. When the distance from the walls corresponds 
to an eddy scale of the radial gradient, a transition takes 
place to a scale independent eddy size turbulent torque 
connecting the two regions at inner and outer boundaries. 
The mean velocity distribution in the above analysis can 
give a reduced shear, Q/Qi ex (i?i/i?)^-^^, rather than 
stable Couette flow. 
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FIG. 3. A schematic drawing of the poloidal magnetic flux 
produced by two coils (left). Superimposed (right) are the flux 
lines from the Maxwell calculation [14 , with the iron shield 
and steel shaft included. This radial flux crosses the high 
shear of the Couette flow producing the enhanced toroidal 
field. The Br field from the Maxwell calculations agrees with 
the calibrated probe measurements to 10%. 



Fig. 3 shows the magnetic field configuration. The 
magnetic field Hall detector probe, internal to the liq- 
uid sodium at the mid-plane, is the primary diagnostic 



of the experiment. It consists of 6 multiple, 3-axis, mag- 
netic field Hall effect detectors at the mid-plane in the 
annular space between the two cylinders and contained 
in an aero-dynamically shaped housing. (The fluid fric- 
tion drag produced by this housing, primarily Ekman 
flow, is estimated to be ^ 0.1 of the end- wall Ekman 
torque.) The l]-gain was then measured using an ap- 
plied calibrated Br magnetic field as a function of the 
coil currents. Because of the high gain in 5^ and the 
lack of perfect orthogonality of the Hall detectors, the 
measured Br would be expected to be contaminated by 
a small fraction of the much larger 5^ field. 
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FIG. 4. (left) The lower curve shows the applied radial field 
Br ^ 12 G and the upper curves show the measured toroidal 
magnetic field for four over- laid experiments at = 68 Hz. 
The right panel shows the ratio Bcf^/Br. 
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FIG. 5. The top trace shows the toroidal field ^ 750G 
vs. time from an applied radial field of ^ 250 G. The Q-gain 
is now reduced to x3. The bottom traces are the AC motor 
torque and current respectively, each showing a 50% increase 
due to back reaction. The delay of several seconds in the 
back reaction torque corresponds to the spin-down time for 
the Couette flow to reach a new, modified velocity profile. 

Fig 4 confirms that the measurements are repeatable 
by showing four experiments over-laid where the Q-gmn 
ratio of x8 is repeated all with a low bias field of Br — 12 
G. The Bcjy produced by the Couette flow shear is about 
c:^ 8 X Br. Note that the ideal ^l-gain could be as high 
as ~ Rm/^TT x20 [22], but the reduced shear in the 
sodium (i.e., not exactly Couette profile) partially ac- 
counts for the reduction in the (^-gain. We then mea- 
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sured whether an increase in the velocity shear in the 
Couette volume might increase the 17-gain ratio, as one 
might naively expect. The velocity shear is increased by 
slowing the outer cylinder by using the DC motor as a 
generator or brake. Therefore the ratio was in- 

creased by 20%, but the l]-gain remains constant. This 
suggests the gain expected due to shear is balanced by 
the greater flux diffusion due to the increased turbulence 
of the now unstable Couette flow (see Fig 2). Roughly 
100 electronic samples are recorded and averaged dur- 
ing the several seconds recording cycle of each of four 
magnetic experiments. The time variation between each 
run reflects slight changes in the Couette flow relaxation 
time and hence, the angular velocity distribution. The 
repeatability among these four runs as well as several ear- 
lier runs in the previous six months gives us confidence 
that the conclusion of high l]-gain is valid. 

The expected back-reaction should occur at a larger 
value of Br where the dissipation and torque of the en- 
hanced Bfjy affects the flow field and the added torque 
shows up as additional current and torque in the AC 



drive motor. This effect is shown in Fig 5, where, in ad- 
dition, the expected reduction in the l]-gain ratio, from 
x8 to ~ x3 is observed when the applied poloidal field 
is Br — 250G. The back reaction stress will modify the 
Couette flow profile such as to reduce the l]-gain and 
hence reduce the back reaction stress until a new steady 
state is achieved at reduced 5^. Fig 5 shows the mea- 
sured difference in power, lOkW, due to the back reac- 
tion. We estimate that the specific magnetic field energy 
density, B'^/Stt is dissipated at < 1] > x{Vol) where 
the volume of high B^p is estimated as 1/3 the length. 
Vol (L/3)(7r/2)(i?| - Rj). This results in a power of 
^ 8 kW, roughly the measured 10 kW. 

In conclusion, a large l]-gain in low turbulent shear in 
a conducting fluid has been demonstrated. This is likely 
to be the mechanism of the l^-gain of a coherent a — Q 
astrophysical dynamo. Driven plumes are planned to be 
used in the a phase of the NMD experiment. 

This experiment has been funded by NSF, Univ. of 
Calif. & LANL, IGPP LANL, and NMIMT. 
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